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We demonstrate theoretically that a dissipationless spin-current can flow a long distance through a diffu¬ 
sive normal metal by using superconductors interfaced with magnetic insulators. The magnitude of this spin- 
supercurrent is controlled via the magnetization orientation of the magnetic insulators. The spin-supercurrent 
obtained in this way is conserved in the normal metal just like the charge-current and interestingly has a term 
which is independent of the superconducting phase difference. The quantum state of the system can be switched 
between 0 and tt by reversing the insulators from a parallel to antiparallel conflguration with an external held. 
We show that the spin-current is carried through the normal metal by superconducting triplet (odd-frequency) 
correlations and that the superconducting phase difference can be used to enhance these, leaving clear spectro¬ 
scopic flngerprints in the density of states. 


INTRODUCTION 

Using superconductors as active components in spintron- 
ics devices is a research field that has attracted increasing 
activity in recent years m Such a synergy becomes possi¬ 
ble both due to the special behavior of spin-polarized quasi¬ 
particles in superconductors, featuring extremely long spin 
lifetimes and spin relaxation lengths Em, and because su¬ 
perconducting Cooper pairs can become spin-polarized EHa. 
This type of Cooper pairs occur not only in superconductors 
with intrinsic triplet pairing, such as UCe 2 El and its cousins 
URhCe and UCoCe m, but can in fact be artificially engi¬ 
neered in hybrid structures between conventional supercon¬ 
ductors and magnetic materials mini. For samples with 
substantial impurity scattering, which is the experimentally 
most common scenario, these triplet Cooper pairs acquire a 
special property known as odd-frequency symmetry IT^ in 
order to satisfy the Pauli principle. What this means is while 
the Cooper pair wavefunction is symmetric under both an ex¬ 
change of space- and spin-coordinates, it is antisymmetric un¬ 
der an exchange of time-coordinates. This property leads to 
remarkable features such as gapless superconductivity fni . 
anomalous Meissner effects CmUTtI . and the possibility to 
create spin-supercurrents in diffusive ferromagnetic materi¬ 
als (TSl. It is worth mentioning that paramagnetic Meissner 
effects and zero-bias conductance peaks can also originate 
from other types of effects which are not related to uncon¬ 
ventional superconductivity, as shown previously in the con¬ 
text of (i-wave superconductivity in the cuprates |[T9ti2Tl and 
more recently for topological insulators (2^ . The quality of 
the materials and avoiding crashing the STM tip into the sam¬ 
ple are of paramount importance for proper identification of 
unusual types of superconductivity such as the odd-frequency 
one mentioned above. 

In the context of utilizing superconductors for spintronics 
purposes, the possibility of spin-supercurrents in ferromag¬ 
netic materials 1^ has earned the triplet Cooper pairs much 
attention. It is known that in structures featuring inhomo¬ 
geneous magnetic order, such as intrinsically textured ferro- 
magnets like Ho ||24l|25l, or multilayers with several ferro- 
magnets |[26l , triplet supercurrents can arise even when using 


conventional s-wave superconductors which feature spinless 
Cooper pairs. However, it can be difficult practically to con¬ 
trol the magnetization direction of each of the individual lay¬ 
ers when using complex structures as in Ref. 1^ to create the 
spin-supercurrent. A large amount of works have studied how 
triplet supercurrents can arise in various types of structures in¬ 
cluding both weakly and strongly polarized ferromagnets (see 
e.g. Refs. Il77ti39ll ). However, it would be highly desirable to 
create a spin-supercurrent fiowing in a normal (non-magnetic) 
metal with a minimum amount of magnetic elements required 
due to the ensuing simplification in how to control the exis¬ 
tence, and the properties, of the spin-supercurrent. 

Now, it was recently experimentally shown in Ref. 1401 
that by using magnetic insulators (MI) in a superconducting 
spin-valve setup, it was possible to control the superconduct¬ 
ing critical temperature Tc by changing the relative magneti¬ 
zation configuration from parallel (P) to antiparallel (AP) by 
application of an external field, causing an essentially infinite 
magnetoresistance effect. This finding prompted us to pose 
the question: what happens when magnetic insulators are used 
in a Josephson junction that has a normal, non-magnetic metal 


Superconductor 



FIG. 1: (Color online) The proposed setup: a Josephson junction 
with magnetic insulators (Mis) inserted between the superconduc¬ 
tors and the normal metal. The magnetic insulators have a magneti¬ 
zation that due to shape anisotropy is expected to be confined to the 
plane perpendicular to the tunneling direction. The magnetic mo¬ 
ments of the Mis on the left and right side of the junction, yxil and 
niR, may be misaligned and an applied superconducting phase dif¬ 
ference across the junction drives both a charge and spin supercur¬ 
rent. 
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as its weak link? Can this create a spin-polarized dissipation¬ 
less flow and, if so, can such a spin-supercurrent be controlled 
externally? We here demonstrate that the presence of mag¬ 
netic insulators in a Josephson junction offers an interesting 
way to create a conserved spin-supercurrent in a normal metal 
(see Fig. [T]). This spin-flow is controlled with the relative 
misalignment angle between the magnetic insulators. A main 
advantage with this setup, compared to previous proposals us¬ 
ing ferromagnets, is that the magnetic conflguration is easily 
tunable since the magnetization of one single very thin ( 1-2 
nm) magnetic insulator needs to be altered. 

Moreover, the spin-supercurrent is carried by odd- 
frequency Cooper pairs which leave a clear spectroscopic fln- 
gerprint in the density of states. We show that by tuning 
the superconducting phase difference 0, one can qualitatively 
change the nature of the proximity effect from a conven¬ 
tional singlet one to a triplet one in situ. We also show that 
the presence of magnetic insulators not only creates a spin- 
supercurrent, but that it has important consequences for the 
quantum phase of the system which undergoes a dynamic 0 - 
TT transition for a single sample by changing from a P to AP 
conflguration for the insulators. 


THEORY 


We use the quasiclassical theory of superconductivity ED 
l42ll in the diffusive limit, where the physics is described by 
the Green matrix function g of the system which is an 8 x 8 
matrix in Keldysh-Nambu space. It is deflned in terms of the 
retarded, advanced, and Keldysh part of the Green function: 


g = 


gR gK 
0 


In the absence of non-equilibrium effects. 


such as applied voltages and temperature-gradients, it is suf- 
flcient to consider the retarded part = g, which may be 
parametrized conveniently as follows 14^ : 


g = 


^^(1 + 77 ) 
. - 2 V 7 


2 V 7 \ 

-V(i+77)y 


,X = i. 


( 1 ) 


We have defined V = (1 — 77 )“^ for normalization and 
the . 7. operation means complex conjugation and reversal of 
quasiparticle energy. The Ricatti-matrices {7,7} are 2x2 ma¬ 
trices in spin space and the Green function g satisfles the Us- 
adel equation fl4l in the normal metal 


Dd^^igd^g) + i[ep 3 ,g] = 0 . ( 2 ) 


Here, D is the diffusion coefficient of the normal metal, ps = 
diag(l, — 1 ), and 5 is the quasiparticle energy measured rela¬ 
tive the Fermi level. In order to account for the magnetic insu¬ 
lators at the interfaces, we use spin-dependent boundary con¬ 
ditions discussed in Ref. ma. The most important effect of 
the magnetic insulators is the spin-dependent phase-shifts ex¬ 
perienced by quasiparticles scattering at the interface, which 
are described by a parameter to be deflned below. The su¬ 
perconducting regions are described by bulk Green functions 


which for the left and right side of the junction are denoted gL 
and gR, where 

= ^giyS/iX + , 7 = -ig_ys /(1 + lc)e"‘®^', (3) 

with j = {L, i?}. We have introduced 

s = sinhO, c = cosh©, with 0 = atanh(Ao/5), (4) 

where Aq is the magnitude of the superconducting order pa¬ 
rameter. The bulk solution is an excellent approximation for 
low interface transparencies when using large superconduct¬ 
ing reservoirs. We have used the second Pauli matrix and 
the superconducting phase difference across the junction is 
deflned sls 0 = Or — Or- The boundary conditions read: 

‘^dCigdxg = [gL,g] +^G^[ML,g] ata; = 0, 

‘^dCngdxg = [g,gR\ - iG^[MR,g] atx = d, 

(5) 

where (j = Rbj/Rn is the ratio between the normal-state 
barrier resistance on side j and the resistance of the nor¬ 
mal metal, and dcpn/ where is the 

transmission probability for channel n and difn are the spin- 
dependent part of the phase-shifts picked up by particles scat¬ 
tered at the interface. Finally, the matrix Mj describes the 
orientation of the magnetic moment of the magnetic insulator 
on side j, while d is the length of the normal metal. Exper¬ 
imentally, it is likely that the magnetic insulators will have 
exchange flelds lying in the plane perpendicular to the tunnel¬ 
ing direction due to shape anisotropy if one uses a layered 
’pancake’ geometry for the junction. This case, and other 
conflgurations, are covered by us setting the right interface to 
Mr = diag(a^,a^) whereas the left interface is allowed to 
have an arbitrary orientation, i.e. Mr = cos adiag(a^, -1- 
sin0sinadiag(a^, —giy) + cos(/)sinadiag(a^, a^). Here, 0 
is the azimuthal angle in the x^-plane and a is the angle be¬ 
tween the magnetization and the 2 ;-axis. For later use, we de- 
flne the magnetic moments of the insulators on the left and 
right side as and niR. 

The boundary conditions used here can also be extended 
da to include a magnetoresistance term Gmr which ac¬ 
counts for the different transmission probabilities for spin-t 
and spin-| particles. Inclusion of such a term amounts mainly 
to an overall reduction of the superconducting proximity ef¬ 
fect and we have explicitly verifled numerically that the spin- 
supercurrent exists even in its presence. The magnetic mo¬ 
ment associated with each magnetic insulator should be un¬ 
derstood as the net average moment of the interface region, 
since a disordered interface might have an internal magnetic 
structure. Moreover, interfaces in hybrid structures are intrin¬ 
sically accompanied by the lack of inversion symmetry. For 
this reason, spin-orbit effects could arise at the interface and 
modify the spin-dependent scattering at the superconducting 
interface Il46l - l48]| . However, it is currently unknown how to 
incorporate such interfacial spin-orbit scattering in the bound¬ 
ary conditions of quasiclassical theory. Nevertheless, even 
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if such a mechanism existed, the spin-dependent phase-shifts 
due to the magnetic insulators captured by the parameter 
in our work is sufficient to produce triplet Cooper pairs, and 
so we do not expect that a second mechanism that accom¬ 
plishes the same thing (due to spin-orbit interaction) would 
bring about any major changes. 

Finally, we will later on include non-ideal effects such as 
spin-flip scattering due to magnetic impurities and spin-orbit 
impurity scattering to see how they influence the spin- and 
charge-flow as well as the density of states in the system. 
These are accounted for (HJ [42l by adding extra self-energy 
terms in the commutator part of the Usadel equation Eq. 

Magnetic impurities: Sgf = - 

( 6 ) 

Spin-orbit scattering: Ego = — rpsgpsT, 

STso 

Here we have defined the matrix vector r = {rx^Ty^Tz), 
where the components are given by: 

= diag(a^,a*), u = {x,y,z}. (7) 

For future reference, we introduce the normalized strength of 
magnetic impurity and spin-orbit scattering as g^f = l/SAoTsf 
and g^o = I/SAoTso, where Tsf/so are the relaxation times as¬ 
sociated with each type of scattering. 


RESULTS 

Spin-Supercurrent via Magnetic Insulators 

We proceed to discuss how the charge- and spin- 
supercurrents sustained by the system are influenced by the 
presence of the ferromagnetic insulators. In the quasiclassical 
framework, these are given by 

Iq = [ deTr{p3(545)^} (8) 

^ 7—00 

and 

Is = j A£Yv{p:iT^{gd^g)^]. (9) 

Here, Nq is the density of states at the Fermi-level in the 
normal-state, e is the electric charge, h is the reduced Planck 
constant, while A is the interface contact area. For future 
use, we also define the bulk superconducting coherence length 
= ^/D/Aq. In the weak proximity effect regime, we were 
able to find a general analytical result for the supercurrents of 
spin and charge (see Appendix for details). We first briefly 
consider the charge-supercurrent which reads: 


Iq = {Iq ,0 + Iq,i cosaG^G^)sm0, (10) 


where the coefficients Iq. 0 and Iq^i are lengthy expressions 
that depend on junction parameters such as the width d, mis¬ 
alignment angle a, temperature T, and the interface trans¬ 
parencies Cl/r- The charge-current is found to be inde¬ 
pendent of which orientation the magnetic moments have in 
the x^-plane, cj). We see that the presence of magnetic in¬ 
sulators coupled to the superconductors introduces a cos a- 
dependence on the supercurrent, not only tuning its overall 
magnitude, but also changing the quantum ground-state of the 
junction from 0 to tt when 

lQ,icosaG^G^ = -lQ,o. ( 11 ) 

Thus, O-TT transitions can now occur even with a normal metal 
interlayer by changing a, a feature which was also reported 
in the ballistic limit in Ref. EH. To demonstrate that this 
is a robust feature, we have computed the charge supercur¬ 
rent without any assumption of a weak proximity effect, thus 
using the full Riccati parametrization. This is shown in Fig. 
I^d), where the current changes sign at a 0.27r correspond¬ 
ing to the O-TT transition. Further information may be inferred 
from the analytical expression for the charge-supercurrent in 
the Appendix, Eq. ( p^ : as the width d of the junction in¬ 
creases, larger values for the spin-dependent phase-shifts 
are required in order to make the O-tt transition possible. 

Interestingly, there exists not only a superfiow of charge in 
the system, but also of spin. The polarization occurs in the 
direction x uir, and so we find that while = 0, one 
has: 


Is = G^G^ sm(l)sma{Is, 0 + Is,icos9). (12) 

Eq. ( p^ is one of the main results of this work. It is seen that 
the spin-supercurrent vanishes if one only has one magnetic 
insulator, in which case G^ or G^ is zero. Moreover, it is 
proportional to sin a, which shows that it is also absent in the 
P or AP alignment {a = 0, tt). For other angles a, however, it 
is in general present. The coefficients are purely 

real and vanish in the absence of superconductivity (A = 0). 
There exists a simple relation between the components of the 
spin-supercurrent in the x^-plane: 


/| sin (j) 
cos f' 


(13) 


This spin-supercurrent has several remarkable features: first 
of all, it is conserved throughout the normal metal just like the 
charge-current. Secondly, it is long-ranged as it fiows through 
a normal metal without any exchange field. Thirdly, it has one 
component that is independent of the superconducting phase 
difference 0. The other component goes like cos6>, mean¬ 
ing that the total spin-supercurrent satisfies /|(6>) = 

This can be understood physically, since a spin-current is in¬ 
variant under time-reversal symmetry. The latter operation 
transforms 0 {—0) and causes the charge-supercurrent to 

change sign. In the Appendix, we give the full expression of 
the spin-supercurrents including the coefficients and their de¬ 
pendence on the junction parameters. The spin-supercurrent 






4 


vs. misalignment angle a is shown in Fig. [^a)-(c) for differ¬ 
ent junction parameters. 

The spin-supercurrent should be possible to control with a 
weak external field coupling to the magnetic insulators, tuning 
their relative orientation. We underline that our structure, un¬ 
like previous works, does not include any ferromagnetic met¬ 
als. In fact, a conceptually similar experimental structure to 
the one that we propose to use in our manuscript has recently 
been demonstrated in Ref. 1401 . There, the authors inves¬ 
tigated a spin-valve structure consisting of a superconductor 
fianked by two magnetic insulators of slightly different thick¬ 
nesses (both of order a few nm). Applying an external mag¬ 
netic field would then control the magnetization orientation 
of the thinner of these layers. This indicates that our results 
are of experimental relevance using currently available tech¬ 
niques. 

In an experimental setting, the normal metal sample may 
well include some degree of magnetic impurities or spin-orbit 
scattering on impurities. Thus, it is of interest to see how such 
non-ideal effects infiuence the predictions made in this paper. 
Although no tractable analytical expression is accessible in 



ck/tt 

(d) 
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FIG. 2: (Color online) Plot of the spin- and charge-supercurrents in 
the system. We have used = 30 nm and the relative temperature 
T/Tc = 0.02. The interface parameters are set to be equal, — 3 
and C = 2, and the phase difference is the one supporting the critical 
current, d — 7r/2. In (a), we have set d = 20 nm, 0 = 0. In (b), 
we have d = 5 nm and 0 = 7r/4. In (c), we set d = 20 nm and 
0 = 7r/4. As expected, the components of the spin-supercurrents 
are mirror-images of each other in (b) and (c) due to the choice of 
magnetic configuration of the insulators, 0 = 7r/4 [see Eq. 

The charge-supercurrent is independent of 0. As seen, it changes 
sign when going from a = 0 to a = tt, signalling a O-tt transition. 
The normalization constant of the charge-current is do = NoeDAjA 
while for the spin-currents it is do = NohDA/8. The contour plot in 
the bottom panel (d) is the charge-supercurrent in the 0-a plane using 
d = 20 nm, showing the occurrence of the O-tt transition around 
a ~ 0.27r (the dark green region corresponds to the 7r-phase). 


this case, we have computed numerically the charge- and spin- 
supercurrent in the presence of spin-fiip scattering and spin- 
orbit impurity scattering as described by Eq. Interest¬ 
ingly, the dissipationless fiow of charge and spin are affected 
very differently depending on the type of scattering. Consider 
first the charge-supercurrent [top row of Fig. |^. With in¬ 
creasing spin-fiip scattering, the current is monotonically sup¬ 
pressed. However, this is not the case for spin-orbit impurity 
scattering (middle panel). Instead, the O-tt transition point 
vanishes and the current retains its order of magnitude even 
for very large values of g^o- Turning to the spin-supercurrent, 
we find that both magnetic impurity scattering and spin-orbit 
impurity scattering suppress the spin-fiow monotonically. 

The physical origin of the different behavior of the charge- 
and spin-supercurrents when adding magnetic impurities and 
spin-orbit scattering can be traced back to how the singlet 
and triplet superconducting correlations are affected by them 
(531 [54l. It can be demonstrated analytically that the singlet 
component is insensitive to spin-orbit impurity scattering in 
the absence of a magnetic field, as is reasonable since spin- 
orbit scattering respects time-reversal symmetry. On the other 
hand, the triplet component is suppressed as the spin-orbit 
scattering rate increases. Based on this, one can now un¬ 
derstand why the charge-supercurrent evolves different with 
increasing spin-fiip and spin-orbit scattering respectively. In 
the former case, both the singlet and triplet components are 
suppressed, i.e. the total superconducting proximity effect is 
reduced, and the current simply decays montonically. In the 
latter case, however, only the triplet part is suppressed. With 
only the singlet part remaining in the normal metal, there is 
no mechanism to cause a O-tt transition and the sign of the 


(a) (b) (c) 
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FIG. 3: (Color online) Plot of the charge-supercurrent and the com¬ 
ponents of the spin-supercurrent in the presence of spin-flip scatter¬ 
ing due to magnetic impurities and spin-orbit impurity scattering. 
Column (a) corresponds to magnetic impurity scattering (lines corre¬ 
sponding to different values of ^sf), column (b) to spin-orbit scatter¬ 
ing (gso), and column (c) to both present with equal magnitude. The 
parameters are set to d = 5 nm, = 80 nm, G^p = 3, C = 2, 
T/Tc = 0.02, and 6> = 7r/2 







































5 


current remains positive and still of appreciable magnitude. 

The same reasoning can be applied to the spin-current. 
In this case, it is solely the triplet part which is responsible 
for its existence. Since both magnetic impurities and spin- 
orbit scattering suppress the triplet Cooper pairs, the spin- 
supercurrent decays monotonically as the scattering rate in¬ 
creases. It can be seen from Fig. [^that the dependence of the 
spin-supercurrent on the misalignment angle a between the 
magnetic insulator moments goes toward a pure sin a profile 
as ^so becomes large. 

Phase-tunable Triplet Superconductivity 

Besides the appearance of this unusual spin-supercurrent, 
a Josephson junction with magnetic insulator interfaces of¬ 
fers a unique way to control triplet superconductivity as 
we now demonstrate. It has previously been shown that a 
crossover from pure conventional even-frequency pairing to 
odd-frequency pairing is made possible via spin-active in¬ 
terfaces in S|MI|N bilayer junctions ISOl . The pivotal pa¬ 
rameter is the ratio between the spin-dependent phase-shifts 
and the normal-state tunnel conductance, in our notation 
which causes a pure odd-frequency proximity pairing state at 
the Fermi level (s = 0) when > 1 while a pure even- 
frequency state occurs when G^ < 1. Experimentally, this 
is manifested as a large zero-energy peak in the density of 
states when the odd-frequency triplets dominate. Conversely, 
a minigap appears in the spectrum for conventional singlet 
pairing. To observe such an effect, it would be necessary to 
fabricate several samples with a varying ratio G^. One way to 
accomplish this could be to vary the width of the MI interlayer 
in order to tune the tunneling probability. 

Instead, we here show that when spin-active interfaces are 
incorporated in a Josephson junction geometry, the crossover 
from even- to odd-frequency pairing can now be controlled 
by the superconducting phase dijference 0, which in turn is 
determined by the current flowing through the system. This 
offers a new way to induce a triplet proximity effect which 
can be changed in situ within a single sample, simply by vary¬ 
ing 0. The crossover is manifested by the qualitative nature of 
the proximity effect, going from a minigap (conventional sin¬ 
glet proximity effect) to an enhanced low-energy peak (odd- 
frequency triplet proximity effect). In order to probe how the 
change in pairing symmetry is manifested experimentally, we 
here compute the density of states in the normal metal region 
and its phase-dependence numerically which allows us to re¬ 
lax the assumption of a weak proximity effect. The DOS nor¬ 
malized to its normal-state value is obtained from the solution 
of the Ricatti equations via: 

Ar(5,^) =Re{Tr[Ar(l + 77)]}/2. (14) 

To make better contact with experimental measurements, we 
have added a small imaginary part to the quasiparticle ener¬ 
gies, s ^ s + iS where (5/Ao ^ 1, which represents inelastic 
scattering. 


(a) 
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FIG. 4: (Color online) Proximity-induced density of states N{e, 0) 
(normalized to its normal-state value) in the middle of the normal 
metal region. In all cases, we are considering the P configuration 
where both magnetic insulators have moments pointing in the z- 
direction. The strength of the spin-dependent phase-shifts occuring 
at the interfaces are given by (a) G^ — 0.55, (b) G^ — 1.05, and (c) 
G^ — 1.55. We have used the parameters d = 20 nm, — 30 nm, 
and C = 3. 
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We consider the most general case where each supercon¬ 
ducting interface contains a magnetic insulator. The results 
are shown in Fig. where we have focused on the exper¬ 
imentally most accessible configuration with the insulators 
in the P state. Three choices of the strength of the spin- 
active scattering at the interfaces are considered in (a)-(c) with 
= {0.55,1.05,1.55}. It is seen that the nature of the 
superconducting proximity changes qualitatively due to the 
presence of the magnetic insulators. It is known that in the ab¬ 
sence of magnetic elements {G^ = 0), a minigap is induced 
in the normal metal which is largest for 6> = 0 and closes at 
6> = TT. In Fig. Qa), the minigap is prominent at small phase- 
differences 0. However, instead of monotonically closing the 
minigap as 0 is driven towards tt, the density of states becomes 
strongly enhanced at low energies. This feature arises due to 
the odd-frequency symmetry of the triplet Cooper pairs in the 
normal metal ElMIl. When the spin-active scattering, taking 
place at the insulators, becomes stronger in Fig. Qb) and (c), 
the minigap has vanished all-together, leaving behind only a 
clear zero-energy peak in the density of states. 

The qualitative change in the density of states (going from 
fully suppressed to enhanced at low energies) can be seen 
clearly also when keeping the superconducting phase differ¬ 
ence 0 fixed and varying the magnetic configuration a. This is 
shown in Fig. as one changes from a P to AP configuration 
(going from a = 0 to a = tt), the system makes a transition 
from hosting proximity-induced triplet superconductivity to 
singlet superconductivity. Our work thus demonstrates a con¬ 
version between singlet and triplet Cooper pairs in a normal 
metal by tuning either the superconducting phase difference 
or the configuration of the magnetic insulators. This has the 
important advantage that it can be done in situ, as opposed to 
using ferromagnets where e.g. several samples with different 
widths are created to suppress the singlet component relative 
the triplet one. 

It is interesting to note that the density of states for each 
electron-spin is highly non-degenerate and tunable, as shown 
in Fig. This could potentially be utilized in creating large 
thermoelectric effects based on the idea of Ref. f55\ which 
demonstrated that the spin-splitted density of states arising 
in superconductor/ferromagnet hybrids could yield a thermo¬ 
electric figure of merits far exceeding what is obtained in the 
non-superconducting phase. 

Similarly to our treatment of the spin- and charge- 
supercurrent, we also investigate the infiuence of spin-flip and 
spin-orbit scattering on the density of states. In the left col¬ 
umn of Fig. we consider different values for G^ and the 
spin-flip scattering rate. Regardless of the value of G^, in par¬ 
ticular of whether it is smaller than or greater than the critical 
value Gip^c = 1, the infiuence of the superconducting proxim¬ 
ity effect on the DOS is diminished. As discussed previously 
in the context of the charge- and spin-supercurrents, this may 
be understood physically from the fact that magnetic impu¬ 
rities suppress singlet and triplet components alike, such that 
the DOS eventually reverts to its normal-state value for any 
phase-difference and energy. The situation is different when 


considering spin-orbit impurity scattering, shown in the right 
column of Fig. |7] Now, the spectroscopic manifestation of 
the superconducting proximity effect depends on whether we 
are in the singlet-dominated regime G^ < 1 or the triplet- 
dominated regime G^ > 1. For G^ = 0.55, the presence of 
spin-orbit scattering leaves the minigap intact while suppress¬ 
ing the zero-energy peak that emerges as the superconducting 
phase-difference increases. Hence, the superconducting prox¬ 
imity effect remains clearly visible in the DOS. For G^ > 1 
shown in (d) and (f), however, increasing the spin-orbit scat¬ 
tering rate causes the low-energy enhancement of the DOS 
to be absent since the triplet component is suppressed by this 
type of scattering. When applying even stronger values of 
spin-orbit scattering, a clear minigap appears for all the values 
of G^. It should also be noted that the spectroscopic finger¬ 
prints of the superconducting proximity effect are much more 
sensitive toward the presence of magnetic impurity and spin- 
orbit scattering than the charge- and spin-supercurrents, the 
former being suppressed in magnitude faster compared to the 
current at a given value of ^sf/so- 
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FIG. 5: (Color online) Proximity-induced density of states N{e, 0) 
(normalized to its normal-state value) in the middle of the normal 
metal as a function of quasiparticle energy £ and the misalignment 
angle a of the magnetic insulators. In (a), the superconducting phase 
bias is set to ^ = 0, corresponding to zero current-flow. In (b), we 
have 0 — 7r/2, corresponding to the critical current-flow. The other 
parameters are set to d = 20 nm, = 30 nm, G,p = 1.05, C = 5. 
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FIG. 6: (Color online) Spin-resolved density of states (normalized to 
its normal-state value) in the £ — ^ plane. In (b), the total density 
of states is shown whereas in (a) and (c) the spin-| and spin-t con¬ 
tributions are shown, respectively. The parameters used are C = 5, 
Gcp = 0.9, d = 10 nm, = 30 nm, and a = 0. 


(a) (b) 

■1 C 
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FIG. 7: (Color online) Plot of the density of states N(s,0) (nor¬ 
malized to its normal-state value) in the energy-superconducting 
phase difference (e-0) plane. The top row illustrates the case where 
= 0.55 and (a) ^sf = 0.05, (b) Qso = 0.05. The middle row has 
G^ = 1.05 and (c) ^sf = 0.10, (d) Qso = 0.10. The bottom row 
shows Gcp = 1.55 and (e) ^sf = 0.15, (f) Qso = 0.15. The remaining 
parameters are set to d = 20 nm, = 30 nm, and ( = 5, in the P 
configuration. 


DISCUSSION AND CONCLUDING REMARKS 


We here discuss in more detail how our work is related 
to previous findings. In the proposal by Houzet and Buzdin 
1^ . a ferromagnetic trilayer was suggested as the minimal 
structure that would be able to generate a long-ranged triplet 
supercurrent. In another work by Grein et al. 1(3^ . strongly 
polarized ferromagnets with two spin-active interfaces were 
considered, thus in some sense being similar to the trilayer 
system of Ref. 1^ with the exception that the spin-bands 
were now assumed to be completely decoupled in the bulk due 
to the large exchange field. In this case, a spin-supercurrent 
was shown to also be generated. Shomali et al EH stud¬ 
ied the spin-current in a Josephson junction with a ferromag¬ 
netic metal bilayer and it was realized that a long-ranged su¬ 
percurrent in ferromagnets could in fact be generated with 
only two ferromagnets 1^ . albeit only as a higher-order ef¬ 
fect. More precisely, there would be a contribution to a long- 
ranged triplet supercurrent from the second Josephson har¬ 
monic sin(26>). This could make experimental detection dif¬ 
ficult, since the magnitude of the second harmonic latter is 
usually much smaller than the first harmonic, and a very spe¬ 
cific fine-tuning of the junction parameters would be required 
to observe the effect. Spin-supercurrents have also been ana¬ 
lyzed in other types of superconducting structures, including 
magnetic textures such as spirals, and also using intrinsically 
triplet bulk superconductors 158 ^ 6 ^ . Very recently 1641 . spin 
supercurrents in junctions composed of multiband supercon¬ 
ductors coexisting with a spin-density wave state was stud¬ 
ied theoretically. Similar dependencies on the superconduct¬ 
ing phase difference and magnetic misalignment between the 
spin-density waves as in our case was shown, even if the sys¬ 
tem under consideration in this work is physically quite dif¬ 
ferent from Ref. (641. 

The spin-supercurrent reported in our work occurs in the 
first harmonic, i.e. it is not a higher-order effect, meaning 
that it is present without any fine-tuning of parameters in or¬ 
der to suppress the first harmonics in favor of higher ones. 
Moreover, it occurs without use of any ferromagnetic met¬ 
als: the spin superflow takes place in a non-magnetic normal 
metal. There are several different choices for magnetic insu¬ 
lators that can be used in the proposed setup shown in Fig. 

Previous experiments considering superconducting hybrid 
structures have utilized magnetic insulators such as EuO (521, 
EuS (40l . and GdN (56l . The particular choice of magnetic 
insulator also depends on how well it can be grown at the in¬ 
terface between the superconductor and the normal metal. We 
speculate that suitable material combinations to construct our 
setup could be Nb and EuO as the superconductor and mag¬ 
netic insulator, or alternatively NbN and GdN. Concerning the 
phase-dependent density of states in the normal metal, experi¬ 
mental techniques are available for measuring this quantity as 
demonstrated in Ref. I57l in a conventional SNS Josephson 
junction. By integrating the junction in a loop geometry, the 
superconducting phase 0 is then tunable via a minute magnetic 
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flux. Using AFM-spectroscopy, a complete mapping of how 
the density of states evolves spatially through the junction as 
a function of 0 is possible. 

We note that very recently, quasiclassical boundary condi¬ 
tions valid for any strength of the barrier polarization were 
derived ll63]| . This opens the possibility to study theoretically 
systems with very strongly spin-polarized magnetic insulators 
and even half-metallic (fully polarized) ferromagnets. Con¬ 
cerning direct experimental detection of the spin-polarization 
of the supercurrent itself, one possibility could be to, once it 
has been created, inject it into a magnetic material with a weak 
magnetic anisotropy and observe the resulting magnetization 
dynamics as a result of spin-transfer. There have also been 
proposals utilizing optical detection of spin transport through 
non-magnetic metals 1651 as well as electrical detection l66l 
via a spin-current induced Hall effect. Further work is needed 
to clarify precisely if and how this could be possible using 


spin-supercurrents in superconducting structures. 

In summary, we have shown that by integrating supercon¬ 
ductors with magnetic insulators, one arrives at a unique way 
to both create and control triplet superconductivity in a well- 
deflned way with the superconducting phase-difference, and 
to also create a conserved and tunable spin-supercurrent flow¬ 
ing through a normal metal. 
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APPENDIX: DETAILED EXPRESSIONS EOR CHARGE- AND SPIN-SUPERCURRENTS 


We here provide comprehensive results for the analytical expressions of the supercurrents of charge and spin supported by the 
system. In the weak proximity effect, one flnds the following completely general expressions: 


lo = 


jx _ 


-‘-s 


— 


NoeDA 

4 

NohDA 


NohDA 


NohDA 


J detanh(^y^4Re|[2/saa;/s - 2ftda;ft - ftdxft - hd^h] - [.T.]|, 

detanh(^y)4Re|[-(/t + fi)dxft - ftdxift + h)] “ [•••]}> 

J detanh(^y)4Im|[(/^ - fi)d^ft - /A(/t - fi)] + 

i 


d£tanh(^y^4Re|[-/ta3;/t + fidxfi] - 


(15) 


Above, the notation (.7.) means changing the sign of energy and complex conjugate and we deflned the inverse temperature 
P = l/{kBT). It is seen that the spinless singlet correlations fs do not contribute to any of the spin-currents. In the special case 
of a normal metal separating the superconductors, one can work further with the above expressions by inserting the solutions 


(16) 


We then get expressions for the supercurrents which is independent of position: 

Iq = NoeDA J d£tanh^^^Re| 2 i/c[(A>|^A>|^ — B^B^) + 2{AtAt — BfBt) — 2{AsAs — BgBs) + (^ 4 ,^ 4 , “ 


JX _ 


NohDA 


^ d£tanh(^)Re{2ifc[(At + A^At - + (if + “ (^t + 

^“d£tanh(^)Re{2fc[-(At - + {B^ - B^B^ + (it - A^At - {B^ - B^Bt]], 


JZ _ 


NohDA j _ b^BO - {AxAj^ - BtSt)]}- 


(17) 


The coefficients B^} for the singlet and each of the triplet components are determined by the boundary conditions. For 
instance, one finds for the charge-supercurrent that 

lQ=NoeDAAiiO J detanh^^^Rejdi/cF”^ sin(/c(i) sinh^ cosa^ |. (18) 
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upon defining the quantity: 


r = (k^tfCn + + 2(G^)^ cos^ cos^{kd) - {k‘^d?C,LQR 


cos a 




The spin-supercurrent is given by: 


TX _ 
— 


NqHDA sin (j) sin aG^G^ 


J d6tanh^^^Re|4i/cr ^ sm{kd) sinh^ 0 + 02 cos 9 ^ |, 


where we have defined the expressions 


( 19 ) 


( 20 ) 


ai = (2{G^fcos^a- {G^)^ + k^d"^{Cl + Cr) + (G^f'j (^k'^d'^CUR - G^G^cosa^ cos^{kd), (21) 

02 = (^k'^d'^Cl + 2{G^fcos^ a - {G^f^ {k'^d^Cn + {G^f'j cos^{kd) + (^k'^d'^CiCR - G^G^cosaj . (22) 
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